Randomly poled crystals as a source of photon pairs 



Jan Pefina Jr.^ and Jiff Svozih'k"'^ 

^ Joint Laboratory of Optics, Palacky University and Institute of Physics of Academy of Sciences of the Czech Republic, 

17. listopadu 50a, 772 07 Olomouc, Czech RepubliQ 

Generation of photon pairs from randomly poled nonlinear crystals is investigated using analyti- 
cally soluble model and numerical calculations. Randomly poled crystals are discovered as sources of 
entangled ultra broad-band signal and idler fields. Their photon-pair generation rates scale linearly 
with the number of domains. Entanglement times as short as several fs can be reached. Comparison 
with chirped periodically-poled structures is given and reveals close similarity. 
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I. INTRODUCTION 

The first experimentally observed nonlinear optical ef- 
fect, second-harmonic generation, was investigated al- 
ready more than forty years ago by Franken |l|. Since 
that, many other nonlinear effects have been revealed and 
understood even at quantum level. Among them, sponta- 
neous parametric down-conversion (SPDC) with its pro- 
duction of photon pairs belongs to the most fascinating. 
The reason is that two photons comprising a photon pair 
generated in one quantum event of this process are mu- 
tually strongly correlated (entangled) as was discovered 
by Hong, Ou and Mandel in eighties [i]. They used m 
their experiments nonlinear bulk crystals that later be- 
come the most common sources of photon pairs [3] . In or- 
der to observe spontaneous parametric down-conversion 
phase-matching conditions of the interacting three opti- 
cal fields have to be fulfilled. Unfortunately, they cannot 
be naturally fulfilled in many highly nonlinear crystals. 
However, Armstrong [4] has arrived with the concept of 
additional periodic modulation of nonlinear susceptibility 
that has been practically developed in periodical poling of 
nonlinear crystals [sl]. In this concept, wave vector of the 
additionally introduced nonlinear modulation is added to 
the natural phase-matching condition and the so-called 
quasi-phase-matching of the interacting optical fields is 
reached this way. Highly nonlinear materials can be ex- 
ploited since then. We note that shortening of a nonlinear 
medium such that phase-matching conditions loose their 
importance is the only alternative way to periodical pol- 
ing. This approach has been applied in photonic-band- 
gap structures in which optical interference is crucial for 
reaching an efficient nonlinear interaction [^-[llj. 

Periodical poling has occurred to be extraordinarily 
useful. It has provided not only compensation for the 
natural phase mismatch. The ability to tailor the prop- 
erties of emitted optical fields has been revealed soon. It 
is based on using ordered nonlinear domains with variable 
lengths (chirped periodical poling). Presence of domains 
of different lengths in an ordered structure allows an ef- 
ficient nonlinear interaction of fields in a broad spectral 
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range. For example, signal and idler fields with spectra 
several hundreds of nm wide can be generated in chirped 
LiNbOs crystals. On the other hand, domains with dif- 
ferent lengths can also be ordered randomly. A bit sur- 
prisingly, the nonlinear interaction can be efficient even 
in this case sometimes referred as stochastic quasi-phase- 
matching (SQPM). Similarly as ordered poled structures 
the randomly poled structures (RPS) allow spectrally 
broadband nonlinear interaction. It is not surprising that 
efficiency of random structures is worse compared to or- 
dered ones. However, they usually put smaller require- 
ments to polarization properties of the interacting optical 
fields as well as orientation of the nonlinear medium [l^ . 
Also fabrication of RPS is much easier because high pre- 
cision is required in production of ordered periodically- 
poled structures. 

The role of SQPM in ID has already been addressed for 
the process of second-harmonic generation |13|-ll7l and 
the process of difference-frequency generation [1^ . More- 
over full domain random structures allowing SQPM for 
transversal second-harmonic generation have been stud- 
ied in 

Here, we focus our attention to the generation of pho- 
ton pairs in randomly poled ID nonlinear crystals 20]. 
It is shown that spectral properties of photon pairs and 
photon-pair generation rates are comparable in RPS and 
chirped periodically-poled structures (CPPS) [211424] . 
This is very promising for easy production of broadband 
and efficient sources of photon pairs. These broadband 
sources are important, e.g., in metrology (quantum opti- 
cal coherence tomography [2^) or quantum- information 
processing [2^, [23] . We note that broadband photon-pair 
sources can also be constructed using zero group- velocity 
dispersion conditions 28[. However, such conditions can 
be met only for certain pump frequencies considering a 
given material. 

The paper is organized as follows. In Sec. II, a general 
theory of SPDC modified to random structures is pre- 
sented. Photon-pair generation rates and intensity spec- 
tra are studied both for random and chirped structures 
in Sec. HI. Sec. IV is devoted to temporal properties of 
the generated photon pairs. Spatial properties of pho- 
ton pairs are addressed in Sec. V. Temperature behavior 
of the quantities characterizing photon pairs is analyzed 
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in Sec. VI. Sec. VII brings the analysis of fabrication er- 
rors. The role of ordering in chirped structures is studied 
in Sec. VIII. Finally, conclusions are drawn in Sec. IX. 



II. SPONTANEOUS PARAMETRIC 
DOWN-CONVERSION IN POLED NONLINEAR 
CRYSTALS 

The process of SPDC in a nonlinear crystal can be con- 
veniently described by the following interaction Hamilto- 
nian Hi„t MM: 



X E^+\z,t)E'fHz,t)E'l '(z,t)-f H.. 



(1) 



L denotes the crystal length. In Eq. ([T]), the positive- 
frequency part of the pump electric-field amplitude is de- 
noted as E^"*"' and eI"'' (e|~'') stands for the negative- 
frequency part of the signal (idler) electric-field ampli- 
tude operator. The z-dependent second-order suscepti- 
bility tensor x'^' describes poling of the nonlinear mate- 
rial. Vacuum permittivity is denoted as a-nd B means 
the transverse area of the optical fields. Symbol H.c. re- 
places the Hermitian-conjugated term. 

Electric-field amplitudes occurring in Eq. ([Ij can be 
conveniently decomposed into harmonic plane waves with 
frequencies uja and wave vectors ka'. 



Ei H^^i)^-^ f duJaE^a ''{i^a)exp{-ikaZ + iujat), 



a = s,i. 



(2) 



The quantum spectral amplitudes E^ \uja) in Eq. ^ 
can then be expressed using photon creation operators 
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(3) 



c is speed of light in vacuum, h reduced Planck constant, 
and Ha stands for index of refraction of field a. 

First-order perturbation solution of Schrodinger equa- 
tion with the initial vacuum state |vac) in the signal and 
idler fields results in the following two-photon quantum 
state Itfj): 

= / duJs / doji^{uJs,oJi)al{uJs)al{u!i)\vsLc). (4) 



The two-photon spectral amplitude $ introduced in 
Eq. (HI is given as follows: 

^{uJs,UJi) = g{uJs,uJi)E^p'^\LL:s + UJi) 

X F{Ak{Lj,,uj,)), (5) 

where g denotes a coupling constant, g{u!s,u!i) = 
iy/ujgUJi I [2c■K^J Us {LOs)ni (wi)] x'^-* (0) , andi?p^-'(iUp) stands 



for the pump-field amplitude spectrum. The stochastic 
function F introduced in Eq. ([5]) describes SQPM and is 
derived in the form [20]: 



F(Afc) = j ^ dz exp(^Afcz) 



(6) 



Symbol Ak describes the natural phase mismatch for the 
interacting fields, Afc = kp — kg — ki. 

In a periodically-poled structure neighbor domains dif- 
fer in signs of x^^'' nonlinearity and function F defined in 
Eq. ([6]) can be recast into the form: 



F(Afc) = ^(-ir 



dzexpiiAkz). (7) 



Symbol denotes the number of domains and n-th do- 
main extends from z — z„_i to z — Zn- Positions z„ of 
domain boundaries are random and can be expressed as 
Zn = Zn-1 + lo + 6ln (ft = 1, ... , Nl, zq = -L) in our 
model using stochastic Gaussian declinations Sin- The 
basic layer length is determined such that quasi-phase- 
matching is reached, i.e. Iq = ir/Ako, Afco = Ak{uj^,ujf), 
and means central frequency of field a. The random 
declinations Sin are mutually independent and can be 
described by the joint Gaussian probability distribution 
P: 



1 



N, 



■exp(-(5L'^B5L). 



(8) 



Covariance matrix B is assumed to be diagonal and its 
nonzero elements equal I/ct^. Stochastic vector (5L is 
composed of declinations 61^, symbol stands for trans- 
position. Characteristic function G of the distribution P 
in Ec^. ([5]) takes the form: 



N 



G{6K) = (exp(i(5K • (5L))av = '^('^'^i 



(9) 



symbol • means scalar product. Vector (5K of parameters 
of the characteristic function G is composed of elements 
Skj. One-dimensional characteristic function G{Sk) in 
Eq. (O equals e:iq){-a'^Sk'^ /A). 

In order to obtain analytical formulas, we integrate the 
expression for function F in Eq. ([7]) domain by domain 
and modify the contributions of the first and last do- 
mains in such a way that the following simple formula is 
reached: 



F(Afc) = |^^(-l)"exp(zAA:z„). 



Ak 



(10) 



n=0 



As a typical structure contains hundreds of domains, in- 
correct inclusion of fields from the first and the last do- 
mains leads to negligible declinations. The formula in 
Eq. ([TUI can be interpreted such that SPDC occurs only 
in domains with positive susceptibility x*^^^ at doubled 
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amplitude and domains with negative susceptibility x^^^ 
play only the role of a 'linear' filler. This interpretation 
elucidates why a pair of domains having one positive and 
one negative signs of susceptibility x^"^^ forms an elemen- 
tary unit for understanding properties of photon pairs. 



III. PHOTON-PAIR GENERATION RATES 
AND INTENSITY SPECTRA 

Photon-pair generation rate as well as intensity spectra 
can be easily derived from mean spectral density of the 
number of generated photon pairs niuj^^uji). The mean 
spectral density n corresponding to the quantum state 
\^) is defined by the formula 

n{uJs,UJi) = {{^a\{uJs)as{us)a\{uJi)ai[uJi)\'ip))s^^, (11) 

where the symbol ()av means stochastic averaging over an 
ensemble of geometric configurations of an RPS. Assum- 
ing the quantum state |^) written in Eq. ((4]) we arrive at 
the formula: 



n(ws,a;j) = |.g(a;s,a;j)n£^^+Hws + Wi)p 
X (|F(AA:(cj„cj,))r) av ■ 



(12) 



Spectrum Ss of, e.g., the signal field and photon-pair gen- 
eration rate N can then be derived using the expressions: 



N= duJs / dUi72{u}s,U}i). 



(13) 
(14) 



The averaged squared modulus of the phase-matching 
function F as determined by the formula in Eq. (jlO[) can 
be written in the form: 



(|F(Afc)| 



H{5k)[l- H{5k)^ 



i_]Hmi 



[1 - H{5k)] 



(15) 



5k{uJsT'^i) — A.k{LJs,LOi) — A/cg. Symbol c.c. replaces the 
complex-conjugated term. Function H{dk) occurring in 
Eq. is defined as: 



H{5k) = exp[iSklo]G{Ako 



Sk), 



G(Afc) = exp - 



(16) 



The averaged squared modulus {\F{Ak)\'^)av of phase- 
matching function given in Eq. (ITSt determines the aver- 
aged spectral density n and behaves as follows. It holds 
that \H\ < 1 and \H\ = 1 for a fully ordered struc- 
ture, li Sk — in a fully ordered structure, H is real 
(H = 1) and the averaged squared phase- matching func- 
tion (|i^(Afc)p)av reaches its maximum value 4:{Nl + 1)'^ . 
Nonzero phase mismatch Sk shifts H into the complex 
plane which results in lower values of the mean value 



(|F(Afc)p)av The larger the Sk, the smaller the mean 
value (|F(A/(;)p)av Inspection of the formula for H in 
Eq. also shows that the larger the value of the ba- 
sic layer length Iq the faster the decrease of mean values 
(|F(Afc)p)av for given Sk. According to the formula in 
Eq. (I16p the larger the standard deviation u of a ran- 
dom structure the smaller the value of \H\. The de- 
crease of values of |i? | results in an increase of the range 
of values of the phase mismatch Sk in which the aver- 
aged squared modulus (|F(Afc)p)av of phase-matching 
function attains non-negligible values [see the formula in 
Eq. (HID]. 

The formula for averaged squared modulus 
(|i^(Afc)p)av of phase-matching function in Eq. ([15]) 
can be substantially simplified under the assumption 
a2(Afco)2A^L/2> 1: 



mm' 



2Nl 



(Afco + Sk)^ 
1 - G(Afco + Sk) 



1 - 2G(Afco + Sk) cos{Skla) + G(Afco + Sk)^ 



(17) 



Increasing values of phase mismatch Sk lead to greater 
values of the denominator in the fraction on the r.h.s. of 
Eq. p7)) that result in the decrease of values of the av- 
eraged squared modulus (|F(Afc)p)av of phase-matching 
function. On the other hand, increasing values of devia- 
tion (T weaken this behavior. 

For comparison, we consider another type of RPS de- 
fined such that Zn = —L + uIq -\- Sin where Sin is a ran- 
dom declination of the n-th boundary. These 'weakly- 
random' structures are more ordered compared to those 
considered earlier because the change of length of an n- 
th domain is compensated by the change in length of 
an {n + l)-th domain. The averaged squared modulus 
(|i^^~''(Afc)p)av of phase-matching function can be de- 
rived in this case as follows: 



(|F-^(Afc)|2)av = 



\G{Ako + Sk)\^' 



Nr 



(Afc)2 

exp(i(5fc/o) 



X Nl — exp{iSklo)- 



- exp(i(5fcZo) 

— exp[iJfcZoA''L] 
1 — exp(i(5fc^o) 



(18) 



Disorder of the boundary positions manifests itself as a 
filter for the averaged squared modulus (|i^^~''(Afc)p)av 
of phase-matching function, as evident from the expres- 
sion in Eq. ([T8)) . This leads to spectral filtering of the 
spectral density n. This behavior is qualitatively dif- 
ferent from that observed in RPS as described by the 
formula in Eq. (jl5l) indicating broadening of the spectral 
density n with increasing values of the deviation cr. 

Spectral broadening is the most interesting feature of 
ordered CPPS that we consider here for comparison. Po- 
sitions of boundaries in these structures are described by 
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the formula z„ = -L+tiIq + C,' {n- NL/2fll, (' = C/Afco, 
and C denotes chirping parameter. The phase-matching 
function F^^^PiAk) then takes the form 22]: 



X exp 



■exp(iAA:7VL/o/2) 



4AfcC' 



erf(/(7Vi/2)) - erf(/(-iVi/2))] 



(19) 



^C'{^ko + Sk)xlo 



5k 



v/C'(Afco + 5fc) 



The error function erf is defined as erf(a;) = 
2/y^ exp{—y'^)dy. Detailed inspection of the formula 
in Eq. (IT^ reveals that the larger the value of chirping 
parameter C the broader the phase-matching function 
F'='^''P(Afc). 

As an example, we consider spectrally degenerate 
(nearly) collinear down-conversion from periodically- 
poled LiNbOa pumped at the wavelength ~ 775 nm 
by a cw laser beam. The signal and idler photons 
thus occur at the fiber-optics communication wavelength 
\s — Ai = 1.55 fj,m. The crystal optical axis is perpendic- 
ular to the fields' propagation direction and is parallel to 
the vertical direction. All fields are vertically polarized 
and so the largest element of the susceptibility tensor 
is used. The natural phase mismatch for this configura- 
tion is compensated by the basic domain length Iq equal 
to 9.51535 /im. A structure composed of Nl — 700 lay- 
ers is roughly 6.5 mm long and typically delivers 4 x 10^ 
photon pairs per 100 mW of pumping in case of weakly 
random positions of boundaries (small values of variance 
a). 

The most striking feature of RPSs is that the photon- 
pair generation rate N increases linearly with the number 
Nl of domains, independently of the standard deviation 
a of the random positions of boundaries [see Fig. [l{a)]. 
Standard deviation a plays the central role in the deter- 
mination of spectral widths ASg and ASi of the signal 
and idler fields. The larger the value of deviation a the 
broader the signal- and idler-field spectra Ss and Si, as 
documented in Fig. [T]d. This behavior is easily under- 
standable because structures generated with larger val- 
ues of a have statistically a broader spatial spectrum of 
the x*^^) (z) modulation which gives more freedom for the 
fulfillment of quasi-phase-matching conditions. It holds 
that the broader the signal- and idler-field spectra Sg and 
Si the smaller the photon-pair generation rate N [com- 
pare Figs. [UJa) and (b)]. It reflects the fact that con- 
structive interference of fields from different domains is 
enhanced in the area outside the central frequencies a;^ 



and 



whereas this interference is weaken in the area 



around the central frequencies. 

The photon-pair generation rate N increases roughly 
linearly with the number of domains also in the 
case of 'weakly-random' structures described by the 
averaged squared modulus (|i^^~''(Afc)|^)av of phase- 
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FIG. 1. (a) Photon-pair generation rate A'^ and (b) signal-field 
spectral width ASs (FWHM) as functions of the number Nl 
of domains for an ensemble of RPSs with standard deviation 
a equal to m (solid curve), 0.1 xlO~®m (solid curve with 
x), 0.5 xlO~®m (solid curve with A), 1 xlO~^m (solid curve 



with o), and 2 xlO 
arbitrary units. 



m (solid curve with o); a.u. stands for 




FIG. 2. Photon-pair generation rate A'' as it depends on the 
number Nl of domains for an ensemble of 'weakly-random' 
structures with standard deviation a equal to m (solid curve 
with *) and 2 xlO~®m (solid curve). 



matching function in Eq. (|18p . as shown in Fig. The 
greater the standard deviation cr the smaller the photon- 
pair generation rate N. As for the signal-field spectral 
width ASs its values do not practically depend on the 
variance a. 

The behavior of photon-pair generation as observed in 
RPSs can also be found in ordered CPPSs. Also here 
photon-pair generation rate N is linearly proportional 
to the number Nl of domains and spectral widths AS'^ 
and ASi increase with increasing chirping parameter C. 
The main result of our analysis is that this similarity is 
both qualitative and quantitative. For any value of the 
chirping parameter C there exists a value of the standard 
deviation a such that spectral widths ASs and SAi of 
the generated signal and idler fields are the same. More- 
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over (and a bit surprisingly) , also photon-pair generation 
rates N are comparable. This is illustrated in Fig. [3] for 
a chirped structure with Nl = 700 domains. Its signal- 
field spectrum Ss is extraordinarily broad (larger that 
I /im) for sufhciently large values of the chirping param- 
eter C [see Fig. [3lja)]. Signal-field spectra 5s of the same 
width can also be generated from RPSs with sufhciently 
large randomness (i.e., having large values of the devia- 
tion a). Values of the standard deviation a correspond- 
ing to the values of chirping parameter ( are plotted in 
Fig. [21b). Photon-pair generation rates N for RPSs and 
CPPSs are drawn for comparison in Fig.[31[c) in this case. 
Whereas CPPSs give better photon-pair generation rates 
TV for larger values of chirping parameter C, RPSs even 
slightly overcome on average CPPSs for smaller values 
of C- Moreover, the signal-field spectra Sg of individ- 
ual realizations may even be broader which results in 
sharper temporal features. On the other hand, these 
spectra are typically composed of many local peaks (see 
Fig. 4). RPSs thus represent an alternative broadband 
and efficient source of photon pairs with properties com- 
parable to CPPSs. We note, that histograms of domain 
lengths corresponding to RPSs are broader compared to 
those characterizing CPPSs. 



Alternatively, RPSs and CPPSs can be compared un- 
der the requirement of equal photon-pair generation rates 
N. Values of the photon-pair generation rate N decrease 
with the increasing values of chirping parameter C [see 
Fig. [2Ic)]. Transformation curve between standard de- 
viation a and chirping parameter C stemming from the 
requirement of equal generation rates N is monotonous 
and is plotted in Fig. [SJa) in the considered case. The 
corresponding signal- field widths ASs plotted in Fig.[S](b) 
reveal that CPPSs provide broader spectra for the most 
of values of chirping parameter However, the difference 
in spectral widths in CPPSs and RPSs is not dramatic. 



The above presented results for RPSs represent an en- 
semble average. In practical applications, properties of 
individual realizations of a given RPS are naturally im- 
portant. In Fig. [51 we show generation rates N and 
signal-field widths ASg for 10000 realizations of the RPS. 
Histograms of the generation rates N and signal-field 
widths ASs plotted in Figs. Mjo) and are close to 
Gaussian distributions, in accordance with the central 
limiting theorem. 



It holds that the larger the deviation a the smaller the 
relative quadratic fiuctuations SN and 6ASs of photon- 
pair generation rates N and signal-field spectral widths 
ASs, respectively [Sx — \J {{Ax)'^)av/ (a;)av, Ax = x — 
(a;)av] (see Fig. [7]). However, we should note that these 
relative fluctuations are quite large and can even ap- 
proach 40 %. 
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FIG. 3. (a) Signal-field spectral width AS, (FWHM) as a 
function of chirping parameter (. (b) Transformation curve 
between the standard deviation a and chirping parameter 
assuming the same spectral widths ASs- (c) Photon- 
pair generation rate for chirped (A'^'^'^"^^, solid curve with *) 
and random (A^, solid curve) structures and their ratio rjv 
(tn ~ A'^/A'^'^"^, dashed curve) as functions of chirping pa- 
rameter C; Nl = 700. 
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(a) (b) 

FIG. 4. Signal-field spectrum Sa for (a) one typical realization 
of RPS (solid curve) and (b) CPPS (solid curve with *) and 
an ensemble of RPSs (solid curve with o). Spectra Ss are 
normalized such that one photon is emitted; a = 2.1 xlO~®m, 
C = 2.5 X 10V~^ Nl = 700. 
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FIG. 5. (a) Transformation curve between the standard de- 
viation a of RPSs and chirping parameter ( assuming the 
same photon-pair generation rates A^. (b) Signal-field spectral 
widths (FWHM) for random {ASs, solid curve) and chirped 
(AS's'"''^, solid curve with *) structures and their ratio rs 
(rg — Ss/Sf^"^, dashed curve) as functions of chirping pa- 
rameter C; Nl = 700. 



IV. TEMPORAL CORRELATIONS, 
ENTANGLEMENT TIME 

There occurs a strong correlation between possible 

detection times of a signal photon and its twin from 
one photon pair because both photons are generated in- 
side the nonlinear medium at one instant. A finite dis- 
tance between the detection times of both photons is 
a consequence of dispersion properties of the nonlinear 
medium through which both photons at different frequen- 
cies propagate before they leave the crystal. Temporal 
correlations of the signal and idler photons can be conve- 
niently described using a two-photon temporal amplitude 
A defined as: 



(20) 



This amplitude A{ta,ti) gives the probability amplitude 
of detecting a signal photon at time tg and an idler pho- 
ton at time U. 

The simplest experimental method for the determina- 
tion of a typical constant characterizing temporal width 
of the two-photon detection window (entanglement time) 
uses a Hong-Ou-Mandel interferometer. In this interfer- 
ometer, the signal and idler fields mutually interfere on 
a beam-splitter and photons leaving the beam-splitter 
at different output ports are subsequently detected in a 
coincidence-count measurement. The coincidence-count 
rate i?„ depends on a mutual time delay r introduced 
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FIG. 6. (a) Photon-pair generation rates A'^ and signal-field 
spectral widths ASs for 10000 realizations of RPSs (each re- 
alization is depicted as a point), (b), (c) Histograms of rates 
N (b) and widths ASs (c); Nr gives the number of samples 
with given properties. Solid curves in (b) and (c) are best-fit 
Gaussian curves; a = 2.1 xlO^V, Nl = 700. 



between the signal and idler photons. It can be shown 
that temporal extension of the interference part in the 
coincidence-count rate i?„ is proportional to entangle- 
ment time under certain conditions. The coincidence- 
count rate i?„ as a function of relative time delay r is 
described by the following formula: 



Rn{T) = l- q{t), 



(21) 
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(T (10 m) 

FIG. 7. Relative quadratic fluctuations SASs of the signal- 
field spectral width (solid curve) and relative quadratic fluc- 
tuations SN of the photon-pair generation rate (solid curve 
with *) as they depend on standard deviation Nl = 700. 



this case and recast into the following form: 



eyip{iujpT) / duis 
Jo 



(26) 



Rf) — 



X \ (Ak{u}s,ujp- u}s),Ak{ujs,ujp~ 



(27) 



where 

P{r) dtj dt2 

Re[{A{h,t2-T)A*{t2,h-T)),^], (22) 
^0 = ^/ dh dt2{\A{hM)?)..- (23) 



go = g{uj°,uj^). Function T introduced in Eqs. (l26l) and 
()27p incorporates phase-matching conditions into the de- 
scription of temporal properties of photon pairs and is 
defined according to the formula: 

J-(Afc, Afc') = (^^(Afc)^^*(Afc'))av; (28) 

phase-matching function F has been introduced in 

Eq. (nnD. 

Considering RPSs with fluctuations of boundaries de- 
scribed by a Gaussian distribution in Eq. ^ function 
in Eq. takes the form: 



Inserting Eqs. ((20)) and (HJ for the two-photon temporal 
amplitude A and quantum state respectively, into 
Eqs. ([22|l and (|23|) we arrive at the expressions: 



Pir) 



4e2c2 62 Ro 



Re 



OO /•OO 

duis / duji uJsU^i 

"'0 



{^{ujs-, '^i)'J'*(^ii ^s))av exp(ia;ir) exp(-ia;sr) 



(24) 



Ro 



■nti 



2 ^oo 



4e^c262 Jo 



dlOs / duJi UJsUJi{\(^{uJs,Uj{)\^)s 



(25) 



For simplicity, we further assume cw-pumping with 
amplitude ^p at frequency w^, i.e. E'^\ujp) = ^pd{LOp — 
u!p). Formulas in Eqs. (|24l) and (l^5t can be simplified in 



J"(Afc, Afc') = 

J"(Afc, Afc') = 
X ( -i7(Afc) 



AfcAF-^(^^'^^') 

X exp[-i(Afc - Ak')NLlai (29) 



1 - H{Dk) 
1 - H{Dk) 

1 - [-g(Afc)]^^ 
l + H{Ak) 

-I- (Afc i — > -Afc') ; 



iJ(Afc) 



H{Ak) + H{Dk) 



(30) 



Dk ~ Afc — Afc'. Function H occurring in Eq. (pO| has 
been defined in Eq. (IT51) . Symbol (Afc < — > —Afc') in 
Eq. (pO)) replaces the term that is obtained by the indi- 
cated exchange applied to the preceded term inside the 
brackets. 

Considering 'weakly-random' structures, the following 
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form of the function T can be derived: 

l-exp[ii:>fc/o(A^L + l)] 



J-*^'(Afc,Afc') = GiDk)- 
G{Ak)G{Ak') 



1 — exp(z£'fc^o) 

1 - expiiAkNLlo) 
1 — cxp(— iA/c'/o) I 1 — exp(— iA/cZp) 

^ 1 - exp{-iDkNLlo) \ 
1 — exp{—iDklo) j 

+ (Afc < — > -Afc') 



(31) 



On the other hand, function J^^^'"'? attains a simple 
form in case of CPPSs: 



J-'='^"'P(AA:, A/c') = F'=''"'P(AA:)F'='"'P*(Afc') 



(32) 



where the formula for i^'='"''P is written in Eq. ([Tl 

Characteristics of temporal correlations (correlation 
time) between the signal and idler fields can also be ob- 
tained from the measurement of sum-frequency inten- 
sity in a nonlinear medium combining the signal and 
idler photons and having a sufficiently high value of 
T^(2)sum nonlinearity. This process allows us to determine 
the temporal correlation function J*^"™ of intensities of 
the signal and idler fields. Intensity P^"^ of the sum- 
frequency field is given along the formula: 



dt 



(vac|£;(+)(t)4(+)(t-T)|V) 



(33) 



where constant 77*^"™ incorporates the value of 
nonlinearity and quantum detection efficiency. 

The general formula in Eq. (1551) can be recast into 
the following form using the expression for function J- 
in Eq. 



dUp\El+HLUp)\^^ 



/"™(t) 



ielc^B^ JO 



duJs^UJsioJp - LUs) / dL0's\J ^'si'^p - ^'s) 



X g{L0s,L0p - LOs)g*{Lo'^,LOp - w^) exp[-i{ujs - ^^'s)t] 
X T {Ak{ujs,ujp - ujs),Ak{uj'^,ujp - lo'^)) . (34) 

When deriving Eqs. and ([M)) we have assumed that 
the nonlinear medium in which sum-frequency generation 
occurs is ideally phase matched for frequencies present in 
the signal and idler fields. 

A detailed analysis of the expression that gives the 
coincidence-count rate Rn in a Hong-Ou-Mandel inter- 
ferometer reveals an important property: the rate i?„ 
does not depend on phase variations along the signal- 
and idler-field spectra in cw regime. This property is 
frequently referred as nonlocal dispersion cancellation 
[31I [32j . It follows that entanglement time At^*-*^ is 
inversely proportional to spectral widths ASs and ASi 




FIG. 8. Entanglement time At^'^^ (FWHM) as a function 
of chirping parameter ( for ensemble of RPSs with standard 
deviations a derived from the curve in Fig. ^jp (solid curve) 
and CPPSs (solid curve with *); Nl = 700. 



of the signal and idler fields despite their complex pro- 
files. We note that the entanglement time At^'~'^ is 
determined by a temporal extension (FWHM) of the 
coincidence-count interference pattern formed by the rate 
Rnij). Entanglement time Ar^*-*^ thus shortens with 
increasing values of the standard deviation a for RPSs. 



The dependence of entanglement time At 



HOM 



the 



deviation a for an ensemble of RPSs composed of 700 
domains is shown in Fig. [8l We can see in Fig. |8] that 
entanglement times At^*-"^ can be as short as several fs 
for sufficiently large values of the deviation a. This indi- 
cates that temporal quantum correlations can be confined 
into an interval characterizing one optical cycle provided 
that spectral phase variations in the signal and idler fields 
are compensated. Entanglement times At^'^^ of CPPSs 
with the same signal-field spectral widths ASg are plot- 
ted in Fig. [8] for comparison that reveals nearly identical 
entanglement times of both types of structures. Typical 
coincidence-count interference patterns given by i?„ for 
photon pairs generated in both types of structures are 
compared in Fig. O They demonstrate a close similar- 
ity of photon-pairs behavior in a Hong-Ou-Mandel inter- 
ferometer. There occur typical oscillations at the shoul- 
ders of the interference dips. Whereas regular oscillations 
characterize CPPSs, irregular oscillations with larger am- 
plitudes occur for individual realizations of RPSs. How- 
ever, widths of interference dips remain practically un- 
changed for different realizations of RPSs. 

Correlation times At''"™ emerging from sum-frequency 
generation are in general longer than entanglement times 



A- 



.HOM 



observed in a Hong-Ou-Mandel interferometer 



because of a strong phase modulation along the wide 
signal- and idler-field spectra Ss and Si (see Fig. [TU|l . 
Correlation times Ar'^"™ can be even an order of mag- 
nitude greater compared to entanglement times At^'^'^ 
for structures with ultra-broadband spectra. However, 
phase modulation along the spectrum can be compen- 
sated to certain extent which gives shorter correlation 
times At''"™. CPPSs have more regular spectral phase 
behavior (as demonstrated in Fig. llOp and quadratic 
phase compensation is usually sufficient to provide wave- 
packets several fs long. As for individual realizations of 
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T (IQ-^s) 

FIG. 9. Coincidence-count rate Rn as it depends on relative 
time delay t for one realization of RPS (solid curve), CPPS 
(solid curve with *), and an ensemble of RPSs (solid curve 
with o); o- = 2.1 X 10"^, C = 2.5 x 10^"^ Nl = 700. 



90 q 




FIG. 10. Phase of the two-photon spectral amplitude 
${u!B,iiJp — U!s) as it depends on normalized signal-field fre- 
quency 2aJs/i^p for one realization of RPS (solid curve) and 
CPPS (solid curve with*); a = 2.1x10"^, C = 2.5xl0''m-^ 
Nl = 700. 



RPSs, quadratic compensation is less efficient because 
of more irregular phase spectral behavior. Despite this 
values of temporal constants typical for chirped struc- 
tures can be approached [see Fig. [TlT a)]. Provided that 
an ideal phase compensation is reached, both types of 
structures give comparable results [see Fig. [TlT b)] and 
are capable to generate photon pairs with wave-packets 
extending over the duration of one optical cycle. Ex- 
perimentally, pulse shapers have been developed for this 
task and their capabilities in the area of photon pairs 
have already been demonstrated [111 . Comparison of the 
results obtained with quadratic and ideal compensations 
reveals that correlation times At''"™ are approx. two 
times larger if we restrict ourselves to quadratic com- 
pensation. Also the value of quadratic chirp that needs 
compensation differs for individual realizations of RPS. 
This requires an adaptive phase compensator. On the 
other hand phase compensation in case of CPPS can be 
reached in a simpler way, e.g., by inserting a peace of 
suitable material of defined length [34, 35]. Despite this 
RPSs are challenging both for basic physical experiments 
as well as metrology applications. 



0.4 -1 




(b) -r (lO-"s) 



FIG. 11. Sum-frequency field intensity J="™ as a function 
of relative time delay r for one realization of RPS (solid 
curve), CPPS (solid curve with *), and an ensemble of RPSs 
(solid curve with o). In (a) quadratic chirp in the signal- 
field amplitude spectrum is compensated for one realization 
of RPS and CPPS; in (b) complete spectral phase compen- 
sation is assumed. The curves are normalized such that 
/^drr™(r) = 1; cr = 2.1 X lO'^'m, C = 2.5 x 10^"^, 
iVL°°= 700. 



V. CORRELATIONS IN THE TRANSVERSE 
PLANE 

In order to describe spatial properties of the signal and 
idler beams (in the transverse plane) a simple general- 
ization of the model presented in Sec. II has to be devel- 
oped. The inclusion of phase-matching conditions also 
in the directions along the x and y axes and assumption 
of spectrally-flat transverse pump-beam profile Ep± (x, y) 
result in the following separable form of a two-photon 
spectral amplitude $ that additionally depends on radial 
{'^s^ and aziniuthal {(fs, (fi) signal- and idler-field 
emission angles (see Fig. [T^ : 

■X. (^xy{i^s,i^i,'&s,^s,'Oi,Vi), (35) 

where function arises from phase-matching condi- 
tions in the z direction and function <^xy originates in 
phase- matching conditions in the transverse xy plane (see 
Fig.[ll|). Funct ion can be derived in analogy with the 
formula in Eq. (5); 

X F(Afc(Ws, Wi, <^s, i^i, '^0): (36) 
where the stochastic function F has been introduced in 
Eq. Phase-matching conditions in the xy plane give 
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FIG. 12. Geometric scheme for the description of spatial prop- 
erties. Direction of the signal- (idler-)field wave vector fcs (ki) 
is given by radial i9s (i^i) and azimuthal ips (<^i) emission an- 
gles. 



the function ^j:y the following form: 



X exp {iAkxX + iAky-y) (37) 

that includes a pump-beam amplitude profile Ep±{x,y) 
in the transverse plane. Assuming normal incidence of 
the pump beam, the cartesian components of nonlinear 
phase mismatch in Eqs. (1551) and ([57)1 can be written as: 

Akx = ks{ujs) sm{-ds) sm{(ps) + ki{uji) sin(i?i) sm{ipi) 
Aky = fcs(ws) sin(t?s) cos{ips) + h{u!i) sin(z?i) cos^tpi) 
Ak^ = kp{u!s + - ks{ujs)coii{'&s) - ki{uji)cos{'di). 

(38) 

We assume a Gaussian pump-beam transverse 
profile in numerical calculations: Ep±(x,y) — 
l/{TrAxpAyp) exp[—{x/ Axp)^ — (y/Ayp)^] and Axp 
(Aj/p) stands for the pump-beam width along the x (y) 
direction. 

We first pay attention to transverse properties of the 
signal beam only. Its spectral density Ss defined as 

Ss{uJs,'&s,fs) = sin{-ds) J duji j d§i sin(?9i) j dipi 

\^[uJs,iO„^„^s.^^,Vi)? (39) 

depends on the signal-field radial {-ds) and azimuthal ((/?s) 
emission angles. As we study photon-pair emission near 
the coUinear geometry, the signal beam (as well as the 
idler beam) has rotational symmetry along the z axis. 
The dependence of spectral density Ss on signal-field ra- 
dial emission angle ds is shown in Fig. [T31 Investigating 
one realization of RPS we observe a typical 'strip-like' 
behavior depicted in Fig. [T^T a). Fixing the value of ra- 
dial emission angle 'd^ spectrum Ss(ws) is composed of 
many peaks occurring at positions specific for the stud- 
ied realization [compare also with Fig. |3^a)]. Each peak 
changes continuously its central frequency as the radial 
emission angle ds moves. We note that this is typical also 
for layered structures that form band-gaps Q . Averaging 
over many realizations of RPSs smoothes this 'strip-like' 
behavior [see Fig. IT^ b')] and leads to that resembling 




(a) 



0.6 0.8 1.0 1.2 

s, (a.u.) 




(b) 
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FIG. 13. Map of signal-field spectral density Ss as it depends 
on signal-field radial emission angle i?3 for (a) one realization 
of RPS, (b) an ensemble of RPSs, and (c) CPPS; ips^Q deg, 
(7 = 2.1 X 10~V, C = 2.5 X 10''m-^ Nl = 700. 



CPPSs [compare Figs. ITSTb) and [131^]. In these cases 
spectral splitting is observed [s!]. This behavior origi- 
nates in phase-matching conditions along the z direction. 

Integration of spectral densities Ss over the signal-field 
frequency Us gives us densities ti^ of photon-pair num- 
bers that are plotted in Fig. [T3] for the structures studied 
in Fig. [131 Whereas the profile of density Usi'ds) is com- 
plex for one realization of RPS, typical shapes with one 
maximum around a nonzero value of 'ds characterize the 
profiles of density ns{'ds) for an ensemble of RPSs and 
CPPS. 
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FIG. 14. Profile of density Ua of signal-field photon num- 
bers as a function of signal-field radial emission angle 
i9s for one realization of RPS (solid curve), CPPS (solid 
curve with *), and an ensemble of RPSs (solid curve with 
o); n™^'' = max^^ [n3(i?s)]. Plane-wave pumping is as- 
sumed; nsi'dsjips) = / dLUsSs{iUs,i^3,'Ps); ifls = deg, a = 
2.1 xlO^'^m, C = 2.5 X lOSi'^ Nl = 700. 

Correlated area gi of an (idler) photon in a pair rep- 
resents spatial analogy to entanglement time and char- 
acterizes correlations of photon twins in the transverse 
plane. By definition, it gives probability of emitting an 
idler photon into radial emission angle "di and azimuthal 
emission angle (pi provided that its signal twin has been 
emitted in a fixed radial emission angle and azimuthal 
emission angle tps, i.e.: 

(i?j ,(pi;'&s,Vs) = sm{'&s ) sin(i?^ ) 

X J dujs J du}^\<^{ujs,uJi,i^s,'Ps,^2,'Pi)\'^- (40) 

Because we mainly pay attention to beams propagat- 
ing in the vicinity of the z axis, we assume that the 
signal photon is emitted along the z axis {'&s — — 
deg). The correlated area as described by function gi 
in Eq. (j40p then has rotational symmetry and its pro- 
files along the radial emission angle 1?^ for CPPS and 
an ensemble of RPSs nearly coincide, as documented in 
Fig. [TSTa). On the other hand, broader profiles are typi- 
cal for individual realizations of RPSs. These individual 
realizations form compact correlated areas without large 
local peaks (compare with Fig. 2] where spectrum Ss for 
one realization of RPS is plotted). The width Ai^i of 
the correlated area along the radial angle di depends in 
general on phase- matching conditions along the z and 'di 
axes. Thus length of the structure, pump-field (tempo- 
ral) spectral width as well as width of the pump-beam 
waist determine together the width A-di (for more details, 
see (36j). For example, focusing the pump beam, values 
of the radial width Ai?i can be varied nearly by one order 
of magnitude [see Fig. [TSTb)]. This behavior can be eas- 
ily explained by the fact that the more the pump beam is 
focused, the wider its spatial spectrum in the transverse 
plane, and so the weaker the phase-matching conditions 
in this plane. 
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FIG. 15. (a) Radial profile gi{'&i) of the correlated area for 
Aj/p = 1 X 10~^m and (b) radial width Ai?i of the correlated 
area as it depends on pump-beam width Ay^ for one real- 
ization of RPS (solid curve), CPPS (solid curve with *), and 
an ensemble of RPSs (solid curve with o); i?i = i9° S-di-, 
^max _ jnax^ J3i(i?i)]. Radially symmetric pump beam is as- 
sumed, i.e. Aa;p = Ayp; = = deg, (fii = 180 deg, 
■d° = Q deg; cr = 2.1 x 10"V, C = 2.5 x 10^"^ Nl = 700. 



VI. THE ROLE OF TEMPERATURE 



We have seen that an ensemble of RPSs and a CPPS 
have similar properties. This close similarity is preserved 
also when studying temperature dependencies [2lj that 
are in general weak. On the other hand, behavior of 
individual realizations of RPSs manifests a stronger tem- 
perature dependence. However, influence of temperature 
varies from realization to realization. Whereas properties 
of the realization of RPS studied above do not consider- 
ably change with temperature (see Fig. [1^] for the signal- 
field spectral width ASg in the temperature range from 
284 to 300 K), other realizations are more prone to the 
change of temperature. This can be conveniently used for 
efficient temperature modifications of properties of pho- 
ton pairs. We note that these effects have their origin in 
temperature dependence of indexes of refraction [2l|. 
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FIG. 16. Signal-field spectral width AS's as it depends on 
temperature T for one realization of RPS (solid curve), CPPS 
(solid curve with *), and an ensemble of RPSs (solid curve 
with o); o- = 2.1 X 10"'^m, C = 2.5 x lO'^m-^ Nl = 700. 
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FIG. 17. Signal-field spectral width ASs as a function of vari- 
ance (Tcr of the fabrication error for one realization of RPS 
(solid curve) and CPPS (solid curve with *). Averaging over 
the fabrication error was done in 1000 randomly chosen posi- 
tions; (J = 2.1 X lO-'^m, C = 2.5 X 10^"^ Nl = 700. 



VII. THE ROLE OF SMALL RANDOM 
(FABRICATION) ERRORS 

In the fabrication process, a small random error nec- 
essarily occurs [13] • This error is sometimes called a 
duty cycle error and, in general, leads to lowering of 
photon-pair emission rates |18i] . Considering spectral 
widths, they are resistant against this error in uniformly 
periodically-poled crystals [33] [see also Eq. (|18p valid for 
'weakly-random' structures]. On the other hand, spectral 
widths are slightly reduced in CPPSs as documented in 
Fig. [m We can see in Fig. [17] that a (large) fabrication 
error with variance o'er = 5 x 10^ ''m results in the reduc- 
tion of signal-field spectral width AS*.; only by approx. 
10 %. Individual realizations of RPSs are much more 
sensitive to the fabrication error. The observed spectral 
changes depend on individual realizations. As an exam- 
ple, the signal-field spectral width ASg of the sample an- 
alyzed above decreases with the increasing variance Ucr 
of the fabrication error. This is natural, because spec- 
trum of this realization is broader compared to the en- 
semble mean value. We note that it holds also here that 
the narrower the signal-field spectrum, the greater the 
photon-pair generation rate and vice versa. 



VIII. THE ROLE OF ORDERING IN CHIRPED 
PERIODICALLY-POLED STRUCTURES 



The benefit of ordering of individual domains by their 
lengths in CPPS can be quantified as follows. We take 
an ordered structure and divide it into segments con- 
taining d domains. We then randomly position these 
segments in a new artificial structure and finally obtain 
mean values of physical quantities after averaging over 
random positions. In the limiting case of d = 1 we have 
a completely random structure similar to those studied 
above. It can be shown that the signal-field spectral 
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FIG. 18. Signal-field spectral width ASs (solid curve with •) 
and photon-pair generation rate (solid curve with A) as 
functions of segment length d. Averaging over 1000 random 
positions of segments was used in calculations. = 2.5 x 



lO^'m-^ Nl = 700. 



width ASs decreases with the decreasing segment length 
d (see Fig. \TE\\ . This is accompanied by an increase of 
photon-pair generation rate N. This behavior refiects 
the fact that spectra of the fields coming from individ- 
ual domains are combined in a more constructive way in 
the central spectral area with the increasing randomness 
(decreasing value of segment length d). The graph in 
Fig. [m also demonstrates that the requirement for the 
same spectral widths ASs of RPSs and CPPS inevitably 
implies that the histogram of domains' lengths for RPSs 
is broader than that obtained for CPPS. 



IX. CONCLUSIONS 

Properties of photon pairs generated in randomly poled 
structures have been found quantitatively similar to those 
characterizing chirped periodically-poled structures. Es- 
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pccially, ultra broadband signal and idler fields can be 
emitted in randomly poled structures. The accompany- 
ing sharp mutual temporal correlations of the signal and 
idler fields can even reach the duration of a single-photon 
cycle (several fs). Photon-pair generation rates depend- 
ing linearly on the number of domains arc; specific to ran- 
dom structures. Stronger temperature dependencies of 
parameters characterizing photon pairs in random struc- 
tures (their individual realizations) compared to those 
found in chirped periodically-poled crystals have been 
observed. In general, application potential of randomly 
poled structures similar to that of chirped periodically- 
poled structures has been revealed. Contrary to chirped 



periodically-poled structures randomly poled structures 
do not require high precision in their fabrication. This is 
a great promise for the use of randomly poled structures. 
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